On the Ordering Instability of Weakly-Interacting Electrons in a Dirty Metal 
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In a dirty metal, electron-electron interactions in the spin-triplet channel lead to singular correc- 
tions to a variety of physical quantities. We show that these singularities herald the emergence of 
ferromagnetism. We calculate the effective action for the magnetic moment of weakly-interacting 
electrons in a dirty metal and show that a state with finite ferromagnetic moment minimizes this ef- 
fective action. The saddle-point approximation is exact in an appropriate large- A'^ limit. We discuss 
the physics of the ferromagnetic state with particular regard to thermal fiuctuations and localization 
effects. 
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I. INTRODUCTION. 

The presence in a solid of randomly-distributed im- 
perfections causes its electrons to propagate diffusively 
rather than ballistically. As a result, when two electrons 
interact, they tarry in each other's vicinity, and become 
more strongly correlatedi. As shown by Finkelstein^i and 
Castellani, et almi^, this effect leads to an instability of 
the diffusive Fermi liquid fixed point (see also refs. 000) 
which manifests itself through divergent corrections to a 
variety of physical quantities. The instability exists for 
arbitrarily weak coupling and for arbitrarily weak disor- 
der. Its simplest incarnation is in an electron system in 
a weak magnetic field B which suppresses quantum in- 
terference effects and relegates them to subleading order 
in the dimensionless resistance^ t (measured in units of 
h/e^) while interactions make a contribution at leading 
order. Since the relevant interaction is in the spin-triplet 
particle-hole channel, the g factor must be tuned to zero 
if B 7^ so that there is no Zeeman splitting, which 
would cut off the divergences. Though this instability 
was identitified nearly 20 years ago, there is, at present, 
no consensus on its outcome. 

The recent discovery of evidence for a metal-insulator 
transition in two-dimensional electron gases (2DEGs) in 
Si-MOSFETsSiia, GaAsii heterostructures, and Si-Ge^^ 
has led to a re-examination of the problem of interact- 
ing electrons in a random potential. The interesting 
phenomena seen in these experiments occur in a regime 
in which interactions are strong (large Ts) and quan- 
tum interference effects are likely to be important be- 
cause the conductivity is on the order of the quantum 
limit ^ /h and B = 0. Though a number of inter- 
esting theoriesi^ii^ii^iiSiiLiS, have been proposed, none 
have provided a complete explanation of the experimen- 
tal data, which isn't surprising, given that the presum- 
ably simpler problem of weakly-interacting electrons in a 
weak random potential with B 7^ 0, g = remains un- 
solved. It is this simpler problem which is the subject of 
this paper 

Though it is unresolved, the instability caused by the 
triplet interaction is, in fact, a crime with only one plau- 
sible suspect: ferromagnetism. Since the diffusive Fermi 



liquid is unstable, the system must settle into some other 
phase at low temperature; since the instability is in the 
spin-triplet channel, it is presumably one in which the 
spins are ordered. The basic interaction in this problem 
is exchange, which is ferromagnetic. There is no avail- 
able mechanism which generates antiferromagnetic inter- 
actions. To see this, note that the instability occurs when 
the conductance is large, while superexchange occurs in 
or near an insulating state (since it arises when the ki- 
netic energy is treated as a perturbation). There would 
be RKKY interactions between localized spins, such as 
magnetic impurities, if we had any, but we don't. In the 
absence of any antiferromagnetic interactions, we rule out 
spin density-waves (SDWs)i2,, spin glass order, and ran- 
dom singlet phases^*^. The latter two furthermore require 
large effective disorder while the instability of interest 
here is manifest for small effective disorder^"'^ . 

It has been suggested that the instability of the triplet 
interaction signals the formation of 'local moments'—. In 
the weak- interaction, weak-disorder limit, however, the 
instability of the diffusive Fermi liquid only develops at 
large length scales, so the resulting magnetic moments 
are anything but 'local'. Local moments might be ex- 
pected when the energy cost for flipping a spin at short 
length scales is very small. This might occur near the 
Stoner transition of the clean systemS^, where the criti- 
cal susceptibility, x ^ ^/T'^ , is of Curie form for a mean- 
field exponent 7 = 1. There is also a substantial local 
moment regime with Curie susceptibility in the large- 
U Hubbard model (because of the weakness of superex- 
change for large U ) and doped semiconductors (because 
of competing ferromagnetic and antiferromagnetic inter- 
actions). However, these scenarios are not relevant to the 
presumably simpler, though still non-trivial, problem of 
the weakly-interacting, weakly-disordered limit. 

This leaves us with one remaining possibility and, in- 
deed, a number of authors have proposed that the sin- 
gularity of the triplet interaction channel leads to ferro- 
magnetism. The basic underlying physics is simply that 
diffusing electrons interact more strongly because they 
move more slowly and, as a result, the interaction energy 
gained by (at least partially) polarizing the system always 
outweighs the single-particle energy cost, in contrast to 
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a clean system, where the interaction energy gain is a 
constant U which must exceed the single-particle energy 
cost l/Np, i.e. the Stoner criterion {Np is the density of 
states at the Fermi energy). Andreev and Kamene^iSSi 
proposed such a mechanism for finite-size systems by 
computing the disorder enhancement of the interaction 
matrix element for a state with two electrons in single- 
particle states separated by small energies. Heuristi- 
cally, one can say that the two wavefunctions are con- 
centrated in the same minima of the potential. Chamon 
and MuccioloM found a ferromagnetic saddle-point solu- 
tion of the one-loop effective action of Finkelstein's the- 
ory. Their calculation generalizes Stoner's by including 
quantum corrections resulting from interactions between 
diffusion modes. Belitz and Kirkpatrick^fi. advanced the 
same hypothesis when they noted that the phase transi- 
tion which they expected to be driven by the growth of 
the triplet interaction is of the same form as a transition 
to ferromagnetism in a disordered system. 

In this paper, we calculate the effective action for the 
magetization in a weakly-interacting electron system in a 
slightly dirty metal. The effective action is calculated to 
lowest order in the dimensionless resistance t. It exhibits 
the physics described above and is clearly unstable to the 
formation of a ferromagnetic moment. In principle, the 
emergence of some other type of order could intervene 
(this effect would have to be non-perturbative in t) be- 
fore ferromagnetism develops, but the arguments given 
above rule out this possibility. Thus, we conclude that 
in the limit of small t, the singularity of the spin-triplet 
interaction channel leads to ferromagnetism. Our results 
substantiate the conclusions of previous author o^^i^'^i^^ . 

The basic philosophy is simple: since the RG equations 
tell us that the interaction in the triplet channel grows 
in importance as the energy scale is lowered, we try to 
satisfy this term in the action first by allowing the mag- 
netization, X, to develop a non-zero expectation value. 
In order to do this, we compute an effective action 5eff [^] 
for the magnetization X taking Finkelstein's 'non-linear 
cr-model' for interacting electrons in a disordered system 
as the point of departure. Since our goal is to calculate an 
effective action S'cff[X] to lowest order in the resistance 
t, Finkelstein's u-model is merely a convenient way of 
organizing the pertinent diagrams. We believe that it is 
intuitively appealing to follow Finkelstein and work di- 
rectly with a theory which retains only diffusion modes 
and drops the superfluous high-energy electronic degrees 
of freedom, but this is a matter of convenience, and it 
is a option which we can take regardless of whether new 
terms are generated at the two-loop or higher level (i.e. 
regardless of whether the non-linear terms in the theory 
are really related by symmetry, as they would be in a 
cr-model). Our derivation emphasizes the similarity be- 
tween, on the one hand, the RG equation for the triplet- 
channel interaction in a dirty metal and the consequent 
development of ferromagnetism with, on the other hand, 
the RG equation for the Cooper-channel interaction in 
a clean Fermi liquid and the consequent development of 



superconductivity. We believe that this illuminates the 
inevitability of our conclusion. 

Once ferromagnetic order has developed, the system 
becomes essentially one of spinless electrons whose diffu- 
sion is decoupled from the spin waves of the system. In 
this regime, the suppression of the resistivity which ac- 
companies the growth of the triplet interaction parameter 
is gone. Thus, our results indicate that a metallic state, 
if it exists in two dimensions, does not owe its existence 
to the growth of the triplet interaction parameter. 

In two dimensions, our results should be interpreted 
as follows. As we lower the temperature, the system de- 
velops local ferromagnetic order as we cross the mean- 
field transition temperature T^^^ . The spin density is 
controlled by classical thermal fluctuations which pre- 
vent order and lead to a correlation length ^ ~ g(const.)/T 
at low-temperatures. In the regime < T < T*^^, we 
have a system in which charge diffuses while the spins are 
gapped. The localization physics of this system is essen- 
tially the same as that of spinless electrons, so the system 
leaves the diffusive regime and eventually becomes local- 
ized. Thus, our small t analysis cannot be used all the 
way to zero temperature. The trend is towards a ferro- 
magnetic insulator, but the onset of insulating behavior 
can thwart the development of ferromagnetism. Recent 
experiments^^ performed near the metal-insulator tran- 
sition indicate an enhancement of the spin susceptibil- 
ity which is cut off by the onset of insulating behavior. 
The metal-insulator transition is outside of the regime 
of validity of our calculation, but this is, nevertheless, 
qualitatively consistent with our picture. 

In 2 -f e dimensions, when the interaction strength is 
larger than a threshhold value of 0(e), there is a ferro- 
magnetic transition at a finite temperature Tc. If the 
resistance is smaller than a critical value tc = 0{e), then 
the resistance will remain finite down to zero temper- 
ature, where the system is in a metallic ferromagnetic 
state. 



II. THE MODEL 

We begin with a system of interacting electrons mov- 
ing in a quenched random potential V{x) in 2D. The 
imaginary-time action is 

S — Sq + "Sdisordcr + ^int (1) 

where 

So + i^disordcr — 

I tpi^^ix) Ue„ -(- + V{x)\ ipnA^) 

(2) 
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and 

Sint — 



T j riV'Jj,a(k + q)V'm,a(k)-(/'/_^(p-q)f/'i+ 



+ r2 V'i,a(k + q)V'm,a(p)V'/,/3(P " q)V'i+m-n,/3 (k) 

For later convenience, we have written the second term 
in momentum space and have spht the interaction into 
nearly forward scattering, in which electrons with mo- 
menta and spins (k, a), (p,/3) are scattered to (k',a), 
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{p',P), where k' « k, p' « p and exchange scatter- 
ing, in which k' f« p, p' f« k. If the microscopic 
electron-electron interaction is V{q), then the associ- 
ated bare couping constants are Fi = V{0), F2 = 
^ {dO / 2tt)V {2k F cos 9/ 2). As Finkelstein noted, these are 
subject to finite Fermi liquid renormalization in the bal- 
listic regime, even before they begin to flow in the dif- 
fusive regime. We assume that the system is in a weak 
magnetic field, so that we can ignore Cooper scattering, 
for which ki « — k2, ks « — k4. 

In order to perform the average over the quenched ran- 
dom potential, we replicate the theory and then average 
the replic ated functio nal integral over the probability dis- 
tribution y(x)V^(x') = ^^^^^ 5('')(x-x'). Now, we have 



n;a L \ / 



n,m,l;a * 



(4) 



r 



where a = 1,2,. . . ,N is a replica index; eventually, we 
must take the replica limit, AT — > 0. Np is the single- 
particle dcnsity-of-states at the Fermi energy. The result- 
ing effective action now contains four quartic terms: one 
disorder term and two interaction terms. Note that the 
disorder term is non-local in time since it only depends 
on two Matsubara frequencies, and it couples different 
replicas. The interaction terms depend on three Matsub- 
ara frequencies, so they are local in time, and they do 
not couple different replicas. 



III. FINKELSTEIN'S EFFECTIVE FIELD 
THEORY 



va 
nm,cx0 



to). It will be useful to think of 
X, Y both as matrices with indices n, m and also as func- 
tions of a single frequency tUn-m', they have only a single 
replica index, however, because interactions do not mix 
different replicas. Note that Xafj <7fja and Y are essen- 
tially the spin and charge density. Integrating out the 
fermions, we can write the partition function as 



/ 



VQVXVY e 



-S^!([Q,X,Y] 



(5) 



We decouple these quartic terms with three Hubbard- 
Stratonovich fields, Y-^ = F«(n 



m), where Sef{[Q, X, Y] is given by 
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Ses[Q,X,Y]^Yl 



-tr In ( ie„ 



2m 



H + —Q + iTiY + T2X 



ttNf 
4t 



tr (g2 



+ lritr(y^) + ir2tr(x2) 



(6) 



In this expression, we use 'tr' to mean the trace over all of the indices which are not explicitly written. Matsubara 
indices arc treated as ordinary matrix indices except that their summations come with factors of T. 
The saddle-point equations are 



ttNfQ 



1 



irV2 +^ + ^Q + iT^Y + T2X 



(7) 



" ^ 2m 



4 



2m ^ 



1 



li+i^Q + iTiY + T2X 



(8) 



y(m) = iT 



2 m 



r 



n+m,m 



(9) 



The momentum integrals on the right-hand-sides of 
© and © are simply iTTNFQn+7n,n, according to iQ. 
Hence, they have the solution 



<9mn,a/3 = SgU (e„) (5„ 

X^^(m) = 

y"(m) = 



(10) 



Note that the right-hand sides of the second and third 
equations of © and vanish for m 7^ because the 
saddle-point Q is diagonal in Matusubara frequency while 
they vanish for to = because of cancellation between 
positive and negative frequencies. As we will see later, 
for r2 sufhciently large, X = becomes unstable and 
there is an X ^ solution of ©• This is the Stoner 
instability. Upon including low-energy fluctuations of Q 
- diffusion modes - an X 7^ solution arises even for r2 
small, as we will see shortly. 



J 



The low-energy fluctuations about this saddle-point 
are transverse fluctuations of Q, which can be 



parametrized by V^'' 



Q 



1/2 



V 



(1 - v^v) 



1/2 



(11) 



The diagonal blocks of Qnm correspond to Matsubara in- 
dices n > 0,m > and 71 < 0,to < while the upper 
right block corresponds to n < 0, to > 0; the lower left 
block, to n > 0,TO < 0. These fluctuations correspond 
to diffusion of charge and spin. The longitudinal fluctu- 
ations of Q as well as fluctuations oi X, Y are gapped at 
the classical level (i.e. tree-level). 

In order to derive an effective field theory for these 
diffusion modes, we shift Q —> Q ~ 2tTi Y + 2rir2 X in 

m 



1 i 
-tr In ( ie„ + — + /i + —Q 
2m 2t 



^^tr (Q2) _ TrTV^ritr {QY) + iTrTVf Fatr {QX) 



-h M 1 27r ( tT^ ) NpTi I Ti tr (F^) + 1 [ 1 - 2^ ( rT^ ) NpT^ | r2 tr ( 



X' 



(12) 



r 



Formally, the factors of t (T^^) in the final two terms 
are infinite as a result of the unrestricted Matsubara 
sum. However, upon integrating out the massive lon- 
gitudinal modes of Q, we see that they are the first 
terms in the series tT^^ — 27r (tT^^^J^ -I- . . . — 
tTY^J (1 + 2tttTJ2^), which is simply l/27r. Thus, we 
replace the factors of r (Tj^m) t>y l/27r in (|12|l . Alterna- 
tively, we could expand the tr ln(. . .) about the diffusive 
saddle-point and keep only terms up to second order in 
X and Y. We would then obtain the same expressions, 
namely H12|l with the factors of t (Tj^m) replaced by 
l/27r. The coefficient of X^ is now 1 — NfT2, which 
becomes negative for T2 > l/Np'- the Stoner instabil- 
ity. Finkelstein actually expanded the tr ln(. . .) only to 
linear order in X and Y . For small Fi, there is no 
difference, but the Stoner instabilty is missed. 



We now expand the tr ln(. . .) about the diffusive 
saddle-point. Thus, we require that Q is a constrained 
field of the form given in . 



5'cff [Q] - t^Nf j d'^x |Dtr(Vg)^ - AiZii (eQ) 

-Fitr (gy)^ 

^ {I + NfTx 



iF2tr {QX) 
Fi 



ttNf 

NfT2) 



■ tr 



ttNf 



tr (X^ 



(13) 



The diffusion constant D is given by D = Vft/2. 
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Following Finkelstein, we integrate out X, Y and ob- 
tain an effective action for the diffusion modes. 

S^s[Q] = t^Nf j (fx |L>tr(VQ)^ - 4itr (eQ) 

(14) 



Note that ordinary matrix multiplication rules reflect the 
non-locality in time of the first term (the 'disorder term') 
in this action. The Fi 2 interaction terms do not in- 
volve matrix multiplication and are, consequently, local 
in time. In this expression, F1.2 — ri_2/(l ± NfTi.2)- 
These corrections to Fi^2 follow from our retention of 
the and terms which Finkelstein drops. For ri^2 
small, Fi 2 = Fi 2, and Finkelstein's effective action is 
recovered. 

Initially, the coefficient of the tr (eQ) term is 1, as in 
but quantum corrections cause it to flow, so we have 
followed Finkelstein in introducing the coupling Z . 



IV. RENORMALIZATION GROUP FLOWS 

The one-loop RG equations associated with this action 
are 



dt 




4 






d£ 


1 2 


dZ 








= t(-Fi 


+ 2F2 


m 






dVi 








= t(r2^ 


4) 


m 




dT2 




2r| 








de 







In 



Z2 
Z 



Zi 



2Fi - T2 



In 



(15) 



where t is the dimensionless resistance, t — 1/{2ttD)\ 
Zi ^ Z - 2Fi + F2; Z2 = Z + F2; and we have set 
■kNp = 1 in order to avoid clutter. 

For sufficiently large interaction strength, the correct 
starting point is presumably a Wigner crystal, which 
corresponds to a radically different saddle-point of the 
effective action Thus, even though these RG equa- 
tions were computed to all orders in Fi^2 for t small, they 
should be understood as holding only for the vicinity of 
the diffusive saddle-point. 

The RG equations for Z, Fi, F2 can be combined to 
give: 



d 

de 



(Z-2Fi+r2) = 



(16) 



Hence, Z — 2Fi + r2 is a constant. For long-ranged 
Coulomb interactions, the Ward identity implies that 
this constant is zero: Z — 2Fi + r2 = 0. More gener- 
ally, we can introduce 71 = Ti/Z, 72 = T2/Z, so that 
71 = ^(1-1-72) — const. /Z. Since Z grows under the RG 



flow l|15() . the flnal term can be dropped and we have 
7i = (1 + 72)72. In the case of spinless electrons, Z does 
not grow under the corresponding RG flow, so we cannot 
make this simplification; we will later discuss the similar 
case of spin-polarized electrons. 

Making the above substitution, we find that the RG 
equations for 72 and Z are: 



drf2_ 

de 

dZ_ 

m 



t{i 

Zi 



72 j 

1 

~2 



■72 



(17) 



Thus, 72 and Z grow as the scale is increased. Eventu- 
ally, this growth causes these equations to break down. 

If we ignore the flow of then the breakdown scale is 
1 

E ~ (1/T)e '(i+T-2) ; in fact, the scale is somewhat lower 
because t flows. 

As noted in ref. 0, the breakdown of these equations 
is reminiscent of the situation in a clean Fermi liquid, 
where the Cooper channel interaction flows according to 
dV/de — —V^. The breakdown of the latter at a scale 
A ~ Ae^^/^ signals the onset of superconductivity or, 
at least, pairing. Unlike in the case of a clean Fermi 
liquid with a four-fermion Cooper channel interaction, Z 
diverges as well in our problem. However, the divergence 
of Z is a by-product of the divergence of 72 ; In the Cooper 
pairing problem, if the interaction which drives pairing 
could also contribute to mass renormalization, we would 
have the same situation as in the Finkelstein problem. 
Indeed, this is what happens when pairing is driven by 
a gauge field. The development of pairing not only cuts 
off the flow of the interaction, but also the flow of the 
effective mass^i. Similarly, we expect that the physics 
which resolves the divergence of 72 will also resolve the 
divergence of Z. 



V. EFFECTIVE ACTION FOR MAGNETIC 
DEGREES OF FREEDOM: ONE-LOOP 

To see how the breakdown of (|15|l leads to ferromag- 
netism, let us go back to the action (|13|) . Instead of inte- 
grating out X and Y at the paramagnetic diffusive saddle 
point, let's integrate out V and Y in order to obtain an 
effective action S'eff[j'^] for the magnetization X. 

The spin polarization X acts as a Zeeman field. A con- 
stant field aligned along the z direction, X^^ = Xa^^p 
(taken to be rephca-diagonal), gives a gap to V|x and 
V|| but leaves V-^^ and unaffected. To see this, 
note that a constant field couples only to the diagonal 
Matusbara indices Qnn- Expanding (1 — y\/t)i/2 
— (1 — yty)i/2 retaining only the leading terms, 
which are quadratic, we find that 



iT2tr{XQ) = iFaXtr 



(18) 



An imaginary gap ±zr2X in Matsubara formalism cor- 
responds to a real gap in real-time. V^i and have 
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propagators: 

(Km;Ti,iT(q)^„m;iT,n(~q)) = 



Z\uJn„,\+Dq^±iT2X ' ' 

When we integrate out V\i and V^i, infrared divergences 
will be cutoff by X, as the single-particle gap cuts off 
divergences in BCS theory. Since V^f and are un- 
affected by X to this order, they simply contribute an 
unimportant constant to S'off[X]. 

We will integrate out V to one loop to obtain S^s [X] ■ 



This is not quite the same as a calculation to lowest or- 
der in t, but the basic physics is illustrated most sim- 
ply in a one- loop calculation, so we will present it first. 
To maintain further simplicity, let us consider the case 
Fi — 0; the general case is a straightforward but more 
complicated extension, to which we return later. In the 
next subsection, we will present a more careful calcula- 
tion which includes the effects of Fi and includes all of 
the terms which contribute to lowest order in t. 

At one loop, we only need to retain terms in up 
to quadratic order in V 



S[V, X] = ^NFY.j |z? tr (vyj Vy+) + Z|c.„,„|yi„^+t/„„,+ + NpT^) ^ tr {X') 

+ zFa [tr {V„„,,+X{n - m)) + c.c] + [tr (v„,,,+V,l+X{n - m)) + c.c.l + {V+ ^ , X ^ -X) \ (20) 



Here, we have written Vj^ = V^i,Vl = V^^, VI ee V^t, and vl = V^^. 

We now shift X in order to remove the bi-linear coupling between X and V . This is the simplest way to proceed 
with a one-loop computation of S'cff[X]. To go beyond one-loop, it is more convenient to shift V , as we will see later. 

Performing the shift X{1) ^ X{1) - iTY.^^^Vm+i,-m;A for ^ > and X{1) ^ X{1) - iTY.n,A^l+i,n; A fo^' ' < 0, where 
A = ±, we have: 

S[V, X]^^NFY.j^Dil (V^t VF) + Z|^„„,|V;L; AVrnn; A ~ TsV^L; A^n'n'; A 5n-m,n'-rn' 

"tr {y^i, aVL, b^\b X{n~ m)) + c.c] +1(1- NpT^) tr {X') 1 (21) 



+ ir2 



A term is also generated as a result of this shift, but 
it doesn't contribute to S'off[Ar] at one loop, so we drop 
it. The V propagator which follows from H21|l is now 

Fa D^^ {u;„„„ q) D^^ (w„„„ g) 6ab (22) 

where 



D^^ (ujnm, q) 



Z\uJ„rn\+Dq^±iT2X 
1 

Z2\t^nm \ + Dq^ ± iT2X 



(23) 



and Z2 is as defined after H15() . 

Integrating out V , we obtain all one- loop diagrams of 
H21|l . with X treated as a classical background field or, 
simply, a tr ln(. . .) expression containing the V propa- 
gators in the X background field. Later, we will re- 
turn to the possibility of other saddle-points, but for 
now, let us assume, in the name of simplicity, that X is 
static, translationally-invariant, and replica-symmetric. 




+ 



FIG. 1: The one-loop diagram with X-treated as a classi- 
cal background field can be viewed as a sum over one-loop 
diagrams with X insertions (the thick lines on the right). 



For such the effective action takes the form: 
S,s[X] ^Y.j\^^^ ^P^^) ^2 tr 

-I- trln^[Z|w„„| + Dq^ ± iV2X] (S„„/ 5mm' 5aa' 5w 

— F2 5n-m-n'+'m' 5aa' Sbb' Sab^ (24) 

This may be rewritten in terms of the propagators Dq^ 
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S,s[X] ^Y.j\^^- ^F^^) r2tr (^2) + trln 



T2 D^-Df-Sa^ 



X- 



(25) 



r 



When we take the product of the two factors on the 
right-hand side of H25|) . the Dq^D^~ term has a free 
summation over the rephca index b which vanishes in 
the rephca Umit N ^ for rephca-symmetric sohitions 
Xa = X. Dropping this term and using the relation 



x± 



T2\L0nm\DQ^D^^ = D^^ , we have 



x± 



S,s[X]^Y. jl^i^-NpT^) Tatr {X^ 



- trln {D^+D^- - D^+D^-) } (26) 



Evaluating the integral, we find 



S,s[X]=pL'' 



-(i-Tv^ra) r^x' 



Z Zo 



T2X 



- 1 



(27) 



where is the area of the system. 

As a result of the logarithm on the second line, the 
minimum of this effective action is at non-zero X for 
arbitrarily small We will discuss this non-trivial so- 
lution and the saddle-point equation (or gap equation) 
which it satisfies in more detail in the next section. At 
this stage, the important lesson is that, at the one-loop 
level, the weakly-interacting diffusive electron gas is al- 
ways unstable to the development of non-zero X, i.e. fer- 
romagnetism. 



VI. FERROMAGNETIC SADDLE-POINT AT 
ONE-LOOP 

Let us examine more closely the saddle-point equation 
for X in the simplest setting, that of section El If we 
vary equation (|25l) with respect to X (before performing 
the momentum integral and Matsubara sum), then we 
obtain the saddle-point equation: 



l-7Vj.r2 = 



7TNFr2Y, J 



^2 ^2 



JjX+ TjX- 



[Z2\L0nm\ + Dq'Y + ^IX^ [^k-nl + Dq^]' + TIX^ 



(28) 



r 



The integrals on the right-hand-side of the second line 
are, for A = 0, precisely the integrals which determine 
the contributions of diagrams (d) and (e) to the RG equa- 
tion for We will further discuss this point below. 

The right-hand-side of the saddle-point equation is log- 
arithmically divergent in the A ^ limit. Hence, there 
is always a solution 



ToX 



1 



-(l-Afj=.r2)(l+72)/(t7l) 



(29) 



We have taken A = l/r as the ultraviolet cutoff and 
72 — 'KNpT2/Z. This is the non-trivial minimum of the 
action to which we referred at the end of section The 
structure of this saddle-point equation is clearly reminis- 



cent of the BCS gap equation. The integrals on the right- 
hand-side are logarithmically divergent - as expected for 
a marginally relevant interaction. A symmetry-breaking 
order parameter must develop in order to cut them off. 

We can submit our saddle-point solution to an im- 
portant check by verifying its consistency with the RG 
equations. The gap to spin-flip excitations, r2A is an 
observable quantity, and should be independent of the 
ultraviolet cutoff, A. Even though the cutoff is taken 
to be a physical scale l/r, we are free to integrate out 
short-wavelength fluctuations in order to lower the cut- 
off while simultaneously renormalizing the couplings Fi, 
r2, t, and Z. This must occur in such a way that the 
saddle-point solution, r2A, remains unchanged. For this 
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reason, we regard the one-loop calculation of sectionlVlas 
incomplete. It is consistent with an RG equation which 
only includes diagrams (d) and (e) of figureElfor the flow 
of Diagrams (a), (b), and (c) arc not incorporated, 
and neither is the flow of Fi, t, Z. 



(a) 




(d) 



(e) 



FIG. 2: The five basic diagrams which contribute to the renor- 
malization of Solid lines are V propagators. Dotted lines 
are bare X propagators. 



VII. BEYOND ONE-LOOP 

Though the preceeding one-loop calculation captures 
the basic physics of the problem, it does not capture all 
of the order t physics. As we have just noted, the effective 
action (|24|) only includes terms corresponding to two of 
the five diagrams in figure|2]- namely, (d) and (e) - which 
contribute to IjlSI) . 

The apparent contradiction stems from the fact that 
the physics of a one-loop RG equation might only appear 
at higher-loops in a saddle-point calculation. This is sur- 
prising at first glance, but can be understood at a tech- 
nical level by noting that a contribution to the effective 
action of the form t^ln^(A/X) is equally important at the 
saddle-point as a contribution of the form t\n{A/X) since 
the extra logarithm compensates for the extra factor of 
t. A simple example of a model in which this occurs is 
given in the appendix. 

The effects of Fi on the flow of t, Z enter in a 
two-loop calculation of the effective action. Consider the 
diagrams of figure|31 These are two-loop diagrams which, 
when combined with the one-loop diagram of the previ- 
ous section, have the effect of shifting in the effective 
action of the previous section the values of t, and Z 
from their bare values. These diagrams can receive con- 
tributions from non-spin-flip diffusion propagators which 
are not cut off by the magnetization, X. They are charge 
and Sz diffusion modes. Since we expect them to be a 
part of the low-energy theory, we can choose not to inte- 
grate them out or to integrate out modes above a cutoff 
of our choosing. The natural choice is to integrate out 
these modes above F2X so that the resulting effective 
field theory for X, Vn contains only excitations with 
energies less than r2X. (We could alternatively imag- 
ine performing the calculation at finite-temperature so 



that T is effectively the cutoff for both spin-flip and non- 
spin-flip propagators. Such a calculation could be used 
to obtain the transition temperature Tc for the putative 
transition.) This effects the following replacements: 



F2 ^ 

t ^ 

Z -, 

where 

/2(z,Zi,Z2) 



t + fVi f2{z, zi, Z2) In 
Z - tVi In 



A 

rTx 



A 



T2X 



(30) 



2zi , , , 2z2 
■Ji(z,zi) - 



Zl - Z2 
1 



Z\ - Z2 



■ fl{z,Z2) 



in(fa) 

Z2 — z \ z J 



(31) 







FIG. 3: Two- loop contributions to ScfffX] which renormalize 
r2, t, and z. The first diagram renormalizes r2 while the 
second renormalizes t and z. The quadratic vertices are Fi 
interaction vertices. 

Diagrams (a)-(c) of figure |21 cannot contribute directly 
to S'off[A] for static X because the external V lines can- 
not carry zero frequency. The Matsubara indices m, n 
of Vmn necessarily satisfy n > 0, m < 0, so their dif- 
ference cannot vanish. For the same reason, these di- 
agrams cannot be linked with (d) and (e), in the way 
demonstrated in the toy model in the appendix. It is 
possible that contributions related to diagrams (a)-(c) 
enter S'off[A] through even higher-loop contributions or, 
perhaps, when fluctuations in X are included. However, 
this touches upon the unresolved issue of the renormaliz- 
ability of the theory at higher-loops. In the absence of a 
clearer understanding of such higher-loop contributions 
to S'off \X\ , we conclude that the one- loop effective action 
(|27|l with the replacement H3U|I is the most systematic 
result which we can obtain at present. We expect that 
going beyond one-loop in V or including the effects of 
fluctuations of X will quantitatively but not qualitatively 
modify our result. 

In section we removed the bi-linear coupling be- 
tween X and y by shifting X. It may seem more natural 
to remove the bi-linear coupling between X and V by per- 
forming the shift Km Vnm - -Dq ^(g, w„m)A'(n - m). 
If we then integrate out V, we will not find a non- 
trivial replica-symmetric saddle-point (though there may 
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be replica-symmetry-breaking saddle-points) , and the re- 
sulting action will not contain any of the physics of the 
one-loop RG equations. In order to recover this physics, 
we must now go to one-loop in S'cff[X]. Going to one- 
loop in ScS [X] , we recover the saddle-point of our initial 
action of section The relevant diagrams are diagrams 
(d) and (e) of figure 13 but the solid lines are shrunk to a 
point when these diagrams are viewed as 1-loop diagrams 
in Seff[^]. Diagrams corresponding to (a), (b), and (c) 
have dotted lines attached to the solid external legs (via 
the bi- linear coupling between X and V); however, these 
vanish in the static limit for the same reason as above. 
Thus, at this level, there is no advantage in proceeding 
this way, and there is a definite disadvantage in that one 
extra step is necessary. However, it is possible that a 



higher-loop calculation is simpler in this approach. 

The one-loop effective action which we computed in 
section O is actually exact in a generalization of the 
model. The most naive large-iV version of the model 
is not particularly simple. If we simply introduce N fla- 
vors of electrons V'a, ct = 1, . . . , iV, then the Q-matrix 
has two flavor indices, Qab, and the large- iV limit allows 
us to drop all non-planar diagrams. However, we do not 
know how to sum all of the planar diagrams. Instead, we 
consider the following generalization: we suppose that 
there are N flavors of diffuson fields V^, a — 1, . . . ,N. 
This generalization does not have a simple interpreta- 
tion in terms of the original electron variables, but it is 
technically simple. The action (|21|l is generalized to: 



S[V,X] =nNFj2j |l?tr {VV^ VV) + ^kn^ll^L; A.l^.™; A. - ^2V^m;A»Vm'n';Aa5n-m,n' 



tr {V„,,, AaVL X{n - m)) + c.c] + y (1 - NpV^) tr (X^) - 1 T^Y. [t^" KvlVp) + c.c] 



(32) 



Note that there is still only a single X field. The explicit 
factor of N in front of the penultimate term and the fac- 
tor of 1/N in front of the final term ensures that every 
term in the action is 0{N). For this action, the diagram 
which we calculated in section is the sole 0{N) dia- 
gram. All others are 0(1) or smaller. Furthermore, our 
effective action is multiplied by a coefficient N. Thus, in 
the iV — > cx) limit, our one-loop saddle-point calculation 
is exact. 

At first glance, it may seem almost contradictory to 
take an iV — > cx) limit while simultaneously taking a 
replica limit n — > 0. However, the flavor index structure 
in the action H32|l is very different from the replica in- 
dex structure. The latter is chosen to ensure that closed 
fermion loops vanish while the latter selects a preferred 
set of the remaining diagrams. Hence, the two limits are 
consistent. 



VIII. PHYSICS OF THE FERROMAGNETIC 
DIFFUSIVE FERMI LIQUID 

Having concluded that the divergence of the triplet 
interaction amplitude leads to ferromagnetism, the ob- 
vious next question is "so then what?" To answer this, 
we must examine the diffusion of electrons in a ferromag- 
netic metal, and the interaction between diffusion modes 
and spin waves. 

In the ferromagnetic state, we can expand X about 
its perfectly-ordered ground state, which we take to be 



along the z-dircction 



1 — mi 



(33) 



We can also perform a gradient expansion of Scff [X] . The 
result is essentially dictated by symmetry: 



.e,X]=/(e,™A™,+i^(Vm.)V0((Vm)' 



(34) 

where i — 1,2 are the x and y components of X. The 
constant K is the spin stiffness of the ferromagnet. 

We now reinstate the coupling of the order parameter 
to the diffusion modes of the system, as in (|13() or l|20|l . 

S[V,m]=nNFY. J (^Dtr {\/V^ \/V)+Z\cj^^\Vl,V„, 
+ iT2Xo [tr {a,V„^Vl^) + c.c] 
+ iT2 [mi{ujn-7n) tr ((TiV;„„) + c.c.' 



+ ^ / (^eijiLUnmimj + K (yrriif^ (35) 

n '' 

where ai are Pauli matrices and we have dropped higher- 
order terms in rrii. A purist may object that the V's 
should not appear in the action for rrii since the dynam- 
ics for the order parameter was obtained by integrating 
out the y's. If we now calculate with this action, we run 
the risk of double-counting the degrees of freedom of the 
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system. In order to avoid this, we must either calculate 
only to tree-level in V (but allow loops in rrii) or intro- 
duce counter-terms which cancel the final line of . In 
the second approach, by computing loops in V, we will 
cancel the rrii counterterms. 

The physics of this model is not completely manifest in 
the above form (|35|l because a uniform rotation of the fer- 
romagnetic order parameter appears to be an excitation 
- since there are non-derivative terms in rrii ~ when, in 
fact, it isn't. This can be made manifest by rotating the 
quantization axis of the diffusing electrons so that it is 
aligned along the local direction of the order parameter. 
This is done by enacting the variable change 



The action is now 



(36) 



s[v,m] = NFj2j^^DtT (vW vy)+z|c.„™|y„t^y„ 

+ iXor2 [ti {(j.VV'') + C.C.] 
Dtr (VV;J„ V(e,jTO,(cj„_0 [<Tj,V^i])) + c.c. 



+ ^ / (eijiuJnmimj+KC^mi 



(37) 



As noted earlier, there is a gap for spin- flip diffusion 
modes, as given by the second line of l|T7|l . Thus, we can 
ignore these modes at low energies and focus on charge 
diffusion and Sz diffusion; Sx,y are carried by spin waves 
nix.y From the structure of the interaction term on the 
third line of l|37|) . we see that m^.y cannot couple directly 
to charge or Sz diffusion alone, as dictated by rotational 
invariance. A charge or Sz diffusion mode, upon emitting 
or absorbing a spin wave, must then become a spin-flip 
diffusion mode. Thus, this coupling can be ignored at 
energies below the gap for spin- flip diffusion modes. Even 
in the absence of such a gap, the coupling would have 
dimension d/2 in d-dimensional space and be strongly 
irrelevant, as expected for interactions with Goldstonc 
modes - which must occur through derivative couplings 
such as that on the third line of 1)3 7|) . 

Thus, the low-energy physics of the ferromagnetic dif- 
fusive metal state is that of ferromagnetic spin waves de- 
coupled from diffusion in a Zeeman field. The latter has 
been well-studied^, and is similar to the spinless elec- 
tron problem. Interactions suppress the conductivity, as 
may be seen from the RG equations: 



dt 
di 



dZ 



Z 



Z-Th-T 



Z + Vh + Ti ^^[Z-Vh+T^^ 



Ti-r 



H 



z 



dVi dVu 



(38) 



Th is the part of the interaction between Sz = diffusion 
modes which includes a contribution from both the direct 



and exchange channels, Th (QtT*5tT + QllQll)'^ its bare 
value is Th = Ti — 

In two dimensions, these RG equations describe a flow 
towards insulating behavior. However, there is a metal- 
insulator transition in 2 + e dimensions through which 
the system can be tuned by varying the bare resistivity 
through a critical value of 0(e). We examine the situa- 
tion in some more detail in the next section. 

The preceeding discussion is strictly true at zero tem- 
perature. At finite temperatures, thermal fluctuations 
will destroy two-dimensional long-range ferromagnetic 
order, but some signature of the zero-temperature fer- 
romagnetic phase persists. If a finite external magnetic 
field is applied, the magnetization is a scaling function 
of temperature and magnetic field. A large-N calculation 
for an SU{N) magnet theory gives the following expres- 
sion for the scaling function for the susceptibilitj>2&: 



Mo 



In 



{q-e-'^/^)/iq 



(39) 



where r — K/T and h = H/T are the rescaled spin stiff- 
ness and magnetic field respectively, g is a root of the 



equation {q - 



-h/2 



q^e-S'T'-, and Mq is the 



zero temperature magnetization. The susceptibility is 
given by Mo(e*'^'' — l)/(87rrr), which diverges exponen- 
tially as T 0. 

In the limit of weak interactions and weak disorder, 
there is a substantial temperature regime in which there 
is ferromagnetic order with a long - but finite - cor- 
relation length and the resistivity is small. If the di- 
mensionless resistivity, t, is taken to be 0.3 at temper- 
ature 5 K, and r^, the average interspace between elec- 
trons measured in the unit of effective Bohr radius, is 
taken to be 0.1, then the susceptibility is large below 
15 mK. However, the resistivity becomes 0(1) at tem- 
peratures <C 10~^° K. While ferromagnetism sets in 
at very low temperatures, insulating behavior is com- 
pletely academic. In the temperature interval between 
these two scales, we expect the susceptibility to be of 
the above form (|39|l . When is larger, the separation 
between these scales decreases, but we should also be 
more circumspect about using this weak-coupling the- 
ory. Of course, if the magnetic field is removed, we pass 
from the unitary ensemble to the orthogonal one, and 
weak-localization effects should drive the system insulat- 
ing much more rapidly. 



IX. PHASE DIAGRAM IN 2 + e DIMENSIONS 

In 2 + e dimensions, the renormalization group equa- 
tions are modified by the appearance of tree-level dimen- 
sions associated with the spatial dimensionality: 



dt 

dl 
dZ 



2ri -T2 \z 



t{-Ti+2T2) 



11 



By varying we can time the system through a 
paramagnet-ferromagnet transition. 

For t, 72 small, the flows are into the fixed line t = 0. 
This fixed line terminates at r2 = r|*°"''''. For r2 = 0, 
there is a metal-insulator transition a.t t = e. In the 
vicinity of this transition, 72 is relevant. If it flows to 
sufficiently large values, then ferromagnetism will de- 
velop. However, if t flows to zero before this happens, 
then the system will remain paramagnetic. Both possi- 
bilities can occur as soon as 72 becomes nonzero, so the 
ferromagnet-paramagnet boundary which begins on the 
t = axis at r2 = r^'°"°'' terminates at f = 6,72 = 0. 
Ferromagnetism, if it occurs, will have a finite onset tem- 
perature. At lower temperatures, the conductivity may 
be metallic or insulating, depending on t. In the for- 
mer case, we have a ferromagnetic metal. In the latter 
case, the system becomes insulating at lower tempera- 
tures. It is possible that the ferromagnetic order is de- 
stroyed by the onset of insulating behavior. This cannot 
be resolved by the methods in this paper; thus, we have 
left the phase boundary of the paramagnetic insulating 
state intentionally vague. (In fact, it may not be param- 
agnetic, but might instead be a spin glass or some other 
non- ferromagnetic phase.) 



PI 






FI 








FM 




PM 







FIG. 4: Phase diagram in 2 -|- e dimensions. PM is the para- 
magnetic metal phase; FM, the ferromagnetic metal phase; 
Fl, the ferromagnetic insulator. PI is a non-ferromagnetic 
insulator: possibly paramagnetic, but perhaps a spin glass. 



the instability of the weakly-interacting electron system 
in a weak random potential signals the development of 
ferromagnetism. Since there is a single instability occur- 
ring in the weak-interaction limit and it occurs while the 
resistivity is small, we are afforded a great deal of control 
over this problem. The instability bespeaks of the rela- 
tive importance of one term in the action compared to the 
others. We need only satisfy this term - which we do by 
allowing for the development of a ferromagnetic moment 
- without worrying about other terms, since none of them 
are relevant. The demonstrable absence of a competing 
instability is the main justification for a simple mean- 
field theory. This may be formalized in a large-iV gen- 
eralization of the problem. Further credence is given by 
the physical arguments of the introductory section which 
rule out any other resolution of the divergence associated 
with the triplet amplitude. This is, in fact, a remark- 
able state of affairs, given that there is no comparable 
controlled theory for itinerant electron ferromagnetism 
in clean metals. In a clean metal, the Stoner instability 
occurs at strong coupling, where there are other com- 
peting instabilities and no unambiguous way of deciding 
between them. 

It is striking how little is known about the phase di- 
agram of the 2DEG. Even at zero magnetic field, the 
phase diagram as a function of rg and disorder strength 
is only known in the small and large limits for van- 
ishing disorder. Intermediate rg is poorly understood, as 
is finite disorder. We believe that we have, in this pa- 
per, extended terra cognita to include the small , weak 
disorder region. 

We emphasize that our result does not shed too much 
light on the metal-insulator transition which is observed 
in strongly-interacting 2DEGs. The development of fer- 
romagnetism does not lead to a metallic state in two 
dimensions. Once ferromagnetism has developed, weak 
localization (which we have neglected in this paper) and 
interaction corrections both drive the system insulating. 
There is no contradiction with experiments, however, 
since the transition seen there is in the strong-interaction, 
high-resistivity {t ^ 1) regime where our result is not ex- 
pected to be valid. On the other hand, our result does 
rule out theories which conjecture that the divergence of 
the triplet interaction leads to a metallic state at zero- 
temperature in 2D. 
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APPENDIX A: CORRESPONDENCE BETWEEN 
RG EQUATIONS AND AN EFFECTIVE ACTION 

When we transform an action by introducing Hubbard- 
Stratonovich fields and integrating out the original vari- 
ables of the action, the result can be an action which 
looks very different from the original one. However, the 
coupling constants of the two actions are related, and the 
corresponding RG equations should be as well. In the 
simplest cases, this is completely transparent. Consider, 
for instance, a system of N flavors of spinless fermions in 
one-dimension. Suppose that their density is commen- 
surate with the underlying crystal lattice and that their 
interaction respects the U (N) symmetry which mixes the 
fermions: 

^ = V'lfl* {dt + d^) + V'Il* {dt - dx) ipaL 

+ j^{^iR^aL+^iL^^aRy (Al) 

Here, a — 1,2,..., TV. tpaW^aL annihilate fermions of 
flavor a at, respectively, the right and left Fermi points. 
This is the U{N) Gross- Neveu model. Let us consider its 
large- limit. 

First, consider the RG equation for g. There is a single 
one-loop diagram, depicted on the left of fig. |S| All other 
diagrams are smaller by powers of 1 /N. The resulting RG 
equation is: 



di 2tt 



(A2) 





FIG. 5: One- loop diagrams contributing to the RG equation 
for g. 

If we introduce a Hubbard-Stratonovich field A which 
decouples the interaction 

+ A (V^l^^^.i + Vl^Va«) + ^ (A3) 

then we can integrate out ip: 

^eff = ^ A2 - TV tr In {{tdt - td.as + Acti) 6ab) (A4) 

where we have combined V'afli V'aL i'^^o ^ two-component 
spinor on which cti^3 act. A cuts off the infrared loga- 
rithmic divergence of the tr In term. 

£eff = ^A^-^NA' (ln(A/A)-l) 
2g iiT 



N 
2^ 



47r 



In (A/A) 



1 

47r 



A^ 



(A5) 



where A is the ultraviolet cutoff. The second line sug- 
gests that 1^ In (A/A) renormalizes 1/g. Thus, the effec- 
tive action manifestly embodies the physics contained in 
the RG equation ljA2p . This is also clear from the close 
kinship between the one-loop diagram which gives the 
trln(. . .) in l|A4p and the one- loop diagram which yields 
the RG equation HA2p . 

However, strictly speaking, the formal relationship 
between the RG equation and the effective action 
(|A5|) follows from the cutoff-independence of physically- 
measurable quantities, A being the simplest. From ljA5|l . 
we see that the saddle-point value of A is 



A = Ae-2V9 



(A6) 



Differentiating both sides with respect to A and requiring 
that dA/dA vanish, we find: 



A 



dg 
'dA 



El 

'2tt 



(A7) 



which is equivalent to the RG equation (|A2p . 

Now, consider a slightly more complicated situation: 

c - Co{^, ^t) + /:o(x, x^) + A„t(x, xt) 



(A8) 



where £o('0jV'^) is the same as the first line of (|A1|) . 
The RG equation for g now reads: 



dg_ 

de 



9_ 
2n 



V 

2^ 



(A9) 



The second term arises from the diagram on the right- 
hand-side of fig. |31 

Now, consider £off (A, x, X^), where A is introduced to 
decouple the g interaction, as before. We can no longer 
calculate it exactly. However, if we assume that v ^ g, 
we can calculate it to lowest non-trivial order in v. The 
lowest order diagram by which it contributes is the second 
one of fig. El We integrate out -0 completely and integrate 
out all X fluctuations above some cutoff A'. We are free 
to choose A' as we like, of course, but the natural choice 
is A' = A so that our effective action is for excitations 
with energies of order A. The resulting effective action 



C,B = ^A^-^NA^ In (A/A) 
Zg 47r 

- — w2iVA2(ln(A/A))^ 



Air 



+ >Co(x,X^)+A„t(x,X^) (AlO) 
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The last term looks strange. It looks like an O (u^) 
renormalization of 1/g, which would imply a contribution 
of order g^v'^ to dg/d£. Also, this term is proportional 
to (A2ln(A/A)) rather than In (A/A) as we might 
have expected from the second diagram of fig. [3 and the 
RG equation IjAQp . 




FIG. 6: The one- and three- loop contributions to the effective 
action for A which encapsulate the same physics as the one- 
loop RG equation. 

The problem is that we are looking at >Cofr(A), which is 
not observable, rather than the vacuum expectation value 
of A, which is. When we compute A, the two difficulties 



noted above resolve each other. Minimizing the effective 
action, we find that, to lowest non-trivial order in v: 

ln(A/A) = — AH) 

or 

A = Ae-^'^/^s+i'Vs) (A12) 

Requiring A- independence of A, we recover the RG equa- 
tion for g - and also the RG equation for v which we 
didn't compute. 

The lesson is that a one-loop contribution to the RG 
equation for a coupling constant might manifest itself 
in an effective action only through a higher-loop (in the 
above example, three-loop) diagram. In the interacting 
electron problem in a random potential, diagrams (a), 
(b), and (c) could enter 6*05 [X] through three-loop dia- 
grams. 



^ B. L. Altshuler and A. G. Aronov, Solid State Comm. 39, 
115 (1979); Sov. Phys. JETP 50, 968 (1979). 

^ A. M. Finkelstein, Zh. Eksp. Teor. Fiz. 84, 168 (1983), 
Sov. Phys. JETP 57, 97 (1983). 

^ C. Castellani, C. Di Castro, P. A. Lee and M. Ma, Phys. 
Rev. B 30 527 (1984). 

* C. Castellani, C. Di Castro, P. A. Lee, M. Ma, S. Sorella, 
and E. Tabet, Phys. Rev. B 30, 1596 (1984). 

5 D. Belitz and T. R. Kirkpatrick, Rev. Mod. Phys. 66, 261 
(1994). 

^ A. M. M Pruis ken, M. A. Baranov, and B. Skoric, 

cond-mat/9712322 9712323. 
'' C.Chamon, A.W.W. Ludwig, and C. Nayak, Phys. Rev. B 

60, 2239 (1998). 

* Note that Finkelstein uses g to denote resistance, as did 
Chamon, Ludwig, and Nayak, in keeping with his notation. 
However, the more widespread usage is for g to denote 
conductance. To avoid confusion, we adopt the notation of 
Castellani, Di Castro, Lee, and Ma: t denotes resistance. 

® S. V. Kravchenko et ai, Phys. Rev. B 51, 7038 (1995); S. 

V. Kravchenko et ai, Phys. Rev. Lett. 77, 4938 (1996); D. 

Simonian, S. V. Kravchenko, M. P. Sarachik, and V. M. 

Pudalov, Phys. Rev. Lett. 79, 2304 (1997). 

D. Popovic, A. B. Fowler, and S. Washburn, Phys. Rev. 

Lett. 79, 1543 (1997). 
" Y. Hanien et al., Phys. Rev. Lett. 80, 1288 (1998); M. Y. 

Simmons et al., ibid, 80, 1292 (1998); S. J. Papadakis and 

M. Shayegan, Phys. Rev. B, June 15 (1998). 
^2 J. Lam et al., Phys. Rev. B 56, R12741 (1997); P. T. Co- 
leridge et al., ibid, 56, R12764 (1997). 

P. Phillips, Y. Wan, I. Martin, S. Knysh, and 

D. Dalidovich, Nature 395, 253 (1998); P. Phillips, 

|cond-mat/9910122 
" Si and C. Varma, Phys. Rev. Lett. 81, 4951 (1998). 

S. Chakravarty, S. Kivelson, C. Nayak, and K. Voelker, 

Phil. Mag. B 79 859 (1999). 
^'^ B. Altshuler and D. Maslov, Phys. Rev. Lett. 82, 145 



(1999). 

" S. Das Sarma and E. H. Hwang, Phys. Rev. Lett. 83, 164 

(1999) ; ibtd 84, 5596 (2000). 

^® C. Chamon, E. R. Mucciolo, A. H. C. Neto, Phys. Rev. B 
64, 245115/1-12 (2001). 

When the bare interaction is weak, disorder does not al- 
low the former, an SDW, and its necessary concomitant, 
a CDW, to form since the energy scale at which the in- 
teraction becomes strong is small compared to the inverse 
of the elastic scattering time. On the other hand, one can 
imagine that strong interactions could lead to the forma- 
tion of an SDW in the ballistic regime, i.e. at length scales 
shorter than the mean-free path. At longer scales, SDW 
order would be disrupted by impurities, but the correla- 
tion length could be long if the impurity density is low. 
In the problem which we consider, however, interactions 
are assumed to be weak in the ballistic regime. They only 
become strong in the diffusive regime, where SDW order 
will not develop. 
^° R. N. Bhatt and P. A. Lee, Phys. Rev. Lett. 48, 344 (1982). 
What we really mean by this is that spin glass order cannot 
be important at this scale. However, local ferromagnetism, 
which develops at the scale at which the triplet interaction 
becomes strong, could be disrupted at longer scales (at 
which the effective disorder strength is large) so that spin 
glass order - e.g. a cluster spin glass - forms at the longest 
length scales. Furthermore, in an insulating state, caused 
by Wigner crystallization and/or Anderson localization, a 
spin glass could form since there will be competition be- 
tween different multi-spin exchange process which lead to 
effective spin-spin interactions of different signs. 

B. N. Narozhny, I. L. Aleiner, and A. I. Larkin, Phys. Rev. 
B 62, 14898 (2000). 

A. V. Andreev and A. Kamenev, Phys. Rev. Lett. 81, 3199 
(1998). 

C. Chamon and E. Mucciolo, Phys. Rev. Lett. 85, 5607 

(2000) . 



14 



D. Belitz and T. R. Kirkpatrick, Phys. Rev. B 53, 14364 
(1996). 

O. Prus, Y. Yaish, M. Reznikov, U. Sivan, and V. Pudalov, 
|cond-mat/0209142| 



N. Bonesteel, I. A. McDonald, and C. Nayak, Phys. Rev. 
Lett. 77, 3009 (1996). 

N.Read and S.Sachdev, Phys. Rev. Lett. 75, 3509 (1995). 



